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Abstract 

I give a survey of different combinatorial forms of alternating-sign ma- 
trices, starting with the original form introduced by Mills, Robbins and 
Rumsey as well as corner-sum matrices, height-function matrices, three- 
colorings, monotone triangles, tetrahedral order ideals, square ice, gasket- 
and-basket tilings and full packings of loops. (This article has been pub- 
lished in a conference edition of the journal Discrete Mathematics and Theo- 
retical Computer Science, entitled "Discrete Models: Combinatorics, Com- 
putation, and Geometry," edited by R. Cori, J. Mazoyer, M. Morvan, and 
R. Mossed, and published in July 2001 in cooperation with le Maison de 
lTnformatique et des Mathematiques Discretes, Paris, France: ISSN 1365- 



8050, http: //dmtcs . lori . f r.) 



1 Introduction 

An alternating-sign matrix of order n is an n-by-n array of 0's, +l's and — l's 
with the property that in each row and each column, the non-zero entries alter- 
nate in sign, beginning and ending with a +1. For example, Figure 1 shows an 
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alternating- sign matrix (ASM for short) of order 4. 
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Figure 1: An alternating- sign matrix of order 4. 

Figure 2 exhibits all seven of the ASMs of order 3. 

+1 \ / +1 \ / +1 \ / +1 \ 
+1 +1 +1 +1 -1 +1 

+1 / \ +1 / \ +1 / \ +1 / 

0+1 \ / +1 \ / +1 \ 
+ 1 o 0+1 0+1 

+1 / \ +1 / \ 0+1/ 

Figure 2: The seven alternating- sign matrices of order 3. 

Matrices satisfying these constraints were first investigated by Mills, Robbins and 
Rumsey [|18|]. The matrices arose from their investigation of Dodgson's scheme 
for computing determinants via "condensation" (see section [H]). The number of 
ASMs of order n, for small values of n, goes like 1,2,7,42,429,7436,..., and it 
was conjectured by Mills et al. that the number of ASMs of order n is given by the 
product 

l!4!7!---(3ra-2)! 
n!(/i + l)!(/i + 2)!---(2n-l)!' 
However, it took over a decade before this conjecture was proved by Zeilberger [ j3"T| ] . 
For more details on this history, see the expository article by Robbins [JZ2"|], the 
survey article by Bressoud and Propp []7|], or the book by Bressoud [§]. 

Here my concern will be not with the alternating- sign matrix conjecture and 
its proof by Zeilberger, but with the inherent interest of alternating- sign matrices 
as combinatorial objects admitting many different representations. I will present 
here a number of different ways of looking at an ASM. Along the way, I will also 
mention a few topics related to ASMs in their various guises, such as weighted 
enumeration formulas and asymptotic shape. Much of what is in this article has 
appeared elsewhere, but I hope that by gathering these topics together in one place, 
I will help raise the level of knowledge and interest of the mathematical commu- 
nity concerning these fascinating combinatorial objects. 
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2 Corner-sum, heights, and colorings 



Given an ASM [a i y)f/=i of order n, we can define a corner-sum matrix (c ; - • 
of order n by putting c i ■ = J\, 



r<] i, J' 



ij=0 

This definition was introduced in |]24|]. 

by-4 



Figure 3 shows the seven corner-sum matrices of order 3 (note that they are 4 
matrices). 
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Figure 3: The seven corner- sum matrices of order 3. 

Corner-sum matrices, viewed as objects in their own right, have a very simple 
description: the first row and first column consist of O's, the last row and last 
column consist of the numbers from to n (in order), and within each row and 
column, each entry is either equal to, or one more than, the preceding entry. 

Note that the seven corner-sum matrices in Figure 3 correspond respectively 
to the seven alternating- sign matrices in Figure 2. I will adhere to this pattern 
throughout, to make it easier for the reader to verify the bijections between the 
different representations. 

Corner-sum matrices can in turn be transformed into a somewhat more sym- 
metrical form. Given a corner-sum matrix (c ; ■ n define h- : = i+ / — 2c, ,. Call 
the result a height-function matrix (see [[□]]). Figure 4 shows the seven height- 
function matrices of order 3 (4-by-4 matrices). 
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Figure 4: The seven height-function matrices of order 3. 

Height-function matrices have a simple intrinsic description: the first row and 
first column consist of the numbers from to n (consecutively), the last row and 
last column consist of the numbers from n to (consecutively), and any two entries 
that are row-adjacent or column- adjacent differ by 1. 

If one reduces a height-function matrix modulo 3, and views the residues 0, 1, 
and 2 as "colors", one obtains a proper 3-coloring of the n + 1-by-rc + 1 square grid 
satisfying specific boundary conditions. Here "proper" means that adjacent sites 
get distinct colors, and the specific boundary conditions are as follows: colors 
increase modulo 3 along the first row and first column and decrease modulo 3 
along the last row and last column, with the color occurring in the upper left. 
Figure 5 shows the seven such 3-colorings of the 4-by-4 grid. 
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Figure 5: The seven colorings associated with the ASMs of order 3. 

Conversely, every proper 3-coloring of that graph that satisfies the boundary 
conditions is associated with a unique height-function matrix [|5|]. 



3 Monotone triangles and order ideals 

Another way to "process" an ASM is to form partial sums of its columns from the 
top toward the bottom, as shown in Figure 6 for a 4-by-4 ASM. In the resulting 
square matrix of partial sums, the zth row has i l's in it and n — i 0's. Hence 
we may form a triangular array whose zth row consists of precisely those values 
j for which the i, j'th entry of the partial-sum matrix is 1 . The result is called a 
monotone triangle [[TSJ] (or Gog triangle in the terminology of Zeilberger [[Hj]). 
Figure 7 shows the seven monotone triangles of order 3. 
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Figure 6: Turning an ASM into a monotone triangle. 
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Figure 7: The seven monotone triangles of order 3. 

One may intrinsically describe a monotone triangle of order n as a triangular 
array with n numbers along each side, where the numbers in the bottom row are 
1 through n in succession, the numbers in each row are strictly increasing from 
left to right, and the numbers along diagonals are weakly increasing from left to 
right. Zeilberger's proof of the ASM conjecture [ f3"T| ] used these Gog triangles and 
a natural generalization, "Gog trapezoids". 

A different geometry comes from looking at the set of ASMs as a distributive 
lattice. Given two height-function matrices (/?■ y07/=o arK * (My)?y=0' we can 
fine new matrices (called the join and meet) whose i, jth entries are max(/z- ., h\ j) 
and min(/z ■ ■ 1 h\ j) , respectively. These new matrices are themselves height-function 
matrices, and the operations of join and meet turn the set of ASMs of order n into 
a distributive lattice L (see [ j26| ] for background on finite posets and lattices). 

The fundamental theorem of finite distributive lattices tells us that L can be 
realized as the lattice of order-ideals of a certain poset P, namely, the poset of 
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join-irreducibles of the lattice L. There is a nice geometric description of the 
ranked poset P. It has ( 1 ) (n — 1 ) elements of rank 0, (2) (n — 2) elements of rank 
1, (3) (n — 3) elements of rank 2, etc., up through (ft — 1)(1) elements of rank 
n — 1 . These elements are arranged in the fashion of a tetrahedron resting on its 
edge. A generic element of P, well inside the interior of the tetrahedron, covers 4 
elements and is covered by 4 elements. 

Using this poset P, we can give a picture of the lattice L that does not require 
a knowledge of poset-theory (also described in p"T|]). Picture a tetrahedron that is 
densely packed with (1) (ft - 1) + (2) (ft - 2) + (3) (ft - 3) + . . . + (ft - 1) ( 1) balls, 
resting on an edge. Carefully remove the two upper faces of the tetrahedron so 
as not disturb the balls. One may now start to remove some of the balls, starting 
from the top, so that removal of a ball does not affect any of the balls below. There 
are many configurations of this kind, ranging from the full packing to the empty 
packing. These configurations are in bijective correspondence with the ASMs of 
order n, and the lattice operations of meet and join correspond to intersection and 
union. 



4 Square ice 

Zeilberger's proof of the ASM conjecture was followed in short order by a simpler 
proof due to Kuperberg [[15]]. Kuperberg's proof made use of a different represen- 
tation of ASMs, the "6-vertex model" of statistical mechanics. This model is also 
called square ice on account of its origin as a two-dimensional surrogate for a 
more realistic (and still intractable) three-dimensional model of ice proposed by 
physicists [|5p. A square ice state is an orientation of the edges of a square grid or 
a finite sub-graph thereof with the property that each vertex other than vertices on 
the boundary has two incoming arrows and two outgoing arrows. Each internal 
vertex must be of one of the six kinds shown in Figure 8 (hence the name "six- 
vertex model"). The markings under the six vertex-types can be ignored for the 
time being. 



+ 1 -1 

Figure 8: The six vertex-types for the square-ice model. 
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As our finite subgraph of the square grid, we will take the "generalized tic- 
tac-toe graph" formed by n horizontal lines and n vertical lines meeting in n 2 
intersections of degree 4, with An vertices of degree 1 at the boundary. We say 
that an ice state on this graph satisfies domain-wall boundary conditions [[14]] if 
all the arrows along the left and right flank point inward and all the arrows along 
the top and bottom point outward. 

Figure 9 shows the possibilities when n = 3. 



Figure 9: The seven square-ice states for n = 3 with domain-wall boundary 

conditions. 

These states of the square-ice model are in bijective correspondence with 
ASMs. To turn a state of the square-ice model on an n-by-n grid with domain- 
wall boundary conditions into an alternating- sign matrix of order n, replace each 
vertex by +1, —1, or according to the marking given in Figure 8. Kuperberg 
was able to give a simplified proof of the ASM conjecture by making use of re- 
sults about the square-ice model in the mathematical physics literature [|14]]. 

An amusing variant of the square-ice model is a tiling model in which the tiles 
are deformed versions of squares that the physicist Joshua Burton has dubbed 
"gaskets" and "baskets", depicted in Figure 10. (To see why the basket deserves 
its name, you might want to rotate the page by 45 degrees, so that the "handle" of 
the basket is pointing up.) 
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Figure 10: A gasket and a basket. 



The gasket and basket correspond respectively to the first and last vertex-types 
shown in Figure 8; the other five vertex-types correspond to tiles obtained by 
rotating the gasket by 90 degrees or by rotating the basket by 90, 180, or 270 
degrees. The directions of the bulges of the four sides of a tile correspond to the 
orientations of the four edges incident with a vertex. Thus, the seven ASMs of 
order 3 correspond to the seven distinct ways of tiling the region shown in Figure 
1 1 (a "supergasket" of order 3) with gaskets and baskets. 



Figure 11: The seven tilings of an order-3 supergasket with gaskets and baskets. 
For another fanciful embodiment of ASMs as tilings, see the cover of [§]. 

5 Symmetric ASMs and partial ASMs 

Some ASMs are more symmetrical than others. More precisely, the eight-element 
dihedral group D 4 acts on ASMs, and for every subgroup G of D 4 there are ASMs 



that are invariant under the action of every element of G. In p3|], Robbins gave 
some conjectures for the number of ASMs of order n that are invariant under par- 
ticular groups G; for most (but not all) of the subgroups G of D 4 , numerical evi- 
dence suggested specific product-formula. Since then, Kuperberg [[TjJ has proved 
some of these, but others remain conjectural. 

At the same time, one may also look at halves (or even quarters or eighths) of 
ASMs — the fundamental regions under the action of the aforementioned groups 
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G — and look at them in their own right, asking, How many partial ASMs are there 
if one limits attention to such a region? There are some interesting phenomena 
here. 

For instance, for Cj,c 2 ,c 3 each equal to +1 or —1, define N(c l ,c 2 ,c 3 ) as the 
number of 4-by-7 partial height-function matrices of the form shown in Figure 12. 

1 2 3 4 5 6 

1 ? ? ? ? ? 5 

2 ? ? ? ? ? 4 

3 3 + Cj 3 3 + c 2 3 3 + c 3 3 

Figure 12: Half of a height-function matrix of order 6. 

Not surprisingly, the eight values of N(c 1 , c 2 , c 3 ) as c l5 c 2 , c 3 vary are not all equal 
to one another. But it is surprising that the four numbers 

#(1,1,1), 

(#(l,l,-l)+#(l,-l,l)+#(-l,l,l))/3, 
(#(l,-l,-l)+#(-l,l,-l)+iV(-l,-l,l))/3, 

and 

#(-1,-1,-1) 

are all equal. More generally, consider (n+ l)-by-(2n+ 1) partial height-function 
matrices of the following form: 
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n + Cj n n + c 2 n n + c 3 

Figure 13: Half of a height-function matrix of order In. 

Here each c- (1 < i < n) is either +1 or — 1 . Let N(c l , c 2 , . . . , c n ) be the number of 
such partial height-function matrices. Then one finds empirically that for every k 
in — n, — (n — 2), — (n — 4), . . . ,n — 4,n — 2,n, the average of N(c l ,c 2 , ■ ■ ■ ,c n ) over 
all vectors (c l5 . . . ,c„) satisfying c 1 + . . . + c n = k depends only on n, not on k. 
Kuperberg has found an algebraic proof of this using the Tsuchiya determinant 
formula [ [271 ] invented for the study of the square ice model, but there ought to be 
a purely combinatorial proof of this simple relation. 
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6 Weighted enumeration 



There are some interesting results in the literature on weighted enumeration of 
ASMs. Here one assigns to each ASM of order n some weight, and tries to com- 
pute the sum of the weights of all the ASMs of order n. A priori it might be 
unclear why this would be interesting, but with certain weighting schemes one 
gets beautiful (and mysterious) formulas, which are their own justification. 

For instance, following [p"g|], one can assign weight x k to every ASM that 
contains exactly k entries equal to — 1 . What is the sum of the weights of the ASMs 
of order n>. When x=l, this is nothing other than ordinary enumeration of ASMs. 
When x = 2, there is a very nice answer [|18jj: the sum of the weights is exactly 
2«(n-i)/2 when x = 3, the answer is more complicated, but it is roughly similar 
in type to the formula for the case x=l, and roughly similar in difficulty; the "3- 
enumeration" formula was first conjectured by Mills, Robbins and Rumsey [ |l^ j 
and was eventually proved by Kuperberg [[R)p (with corrections provided by Robin 
Chapman). No other positive integer x seems to give nice answers. One can also 
assign weight x k to every ASM that contains exactly k entries equal to +1, but 
this is essentially the same weighting scheme, since in any ASM of order n, the 
number of +1 's is always n plus the number of — 1 's. 

More interestingly, one can also use a hybrid weighting scheme in which the 
exponent of x is equal to the number of entries a ■ ■ such that either i + j is even 
and a ; ■ = — 1 or i + j is odd and a- ■ = +1. When x = 2, this too leads to an 
interesting result: the sum of the weights is always a power of two times a power 
of five! [BO] 



One can come up with many open problems by combining the ideas of this 
section and the previous section. Here is one example: Each way of filling in 
Figure 13 (with the c-'s now permitted to vary freely) gives rise to a "half-ASM". 
For instance, the filling 

1 2 3 4 5 6 

1 2 3 4 5 4 5 

2 3 4 3 4 3 4 

3 2 3 2 3 4 3 

of Figure 12 gives rise to the half-ASM 

+1 
+1 
+1 -1 +1 
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which has a single —1. If we assign each rectangular array that arises in this way 

2 

a weight equal to 2 to the power of the number of — 1 's, we get 2 n (Robin Chap- 
man and Theresia Eisenkolbl have independently found nice proofs of this). On 
the other hand, suppose we now permit every entry in the bottom row of Figure 
13 to vary freely (aside from the n on the left and the n on the right). When 
we 2-enumerate half-ASMs of this sort, as a function of n, we get the follow- 
ing sequence of numbers: 2, 20 = 2 2 ■ 5, 896 = 2 7 • 7, 177408 = 2 8 ■ 3 2 • 7 • 11, 
154632192 = 2 15 - 3 - ll 2 - 13, 592344383488 = 2 17 -ll 2 - 13 3 • 17, . . .. Clearly the 
absence of larger prime factors indicates that there is some nice product formula 
governing these numbers. Can someone find the right conjecture? Can someone 
prove it? (For more data of this kind, see: 



http : / / www . math . wise . edu/ ^propp /half -asm 



Late-breaking news: Theresia Eisenkolbl Q10| ] has made progress with the data-set 



and proved a number of theorems about half-ASMs.) 



7 Full packings of loops 

Given an ice state of order n, we can form a subgraph of the underlying tic-tac-toe 
graph by selecting precisely those edges that are oriented so as to point from an 
odd vertex to an even vertex, where we have assigned parities to vertices so that 
each odd vertex has only even vertices as neighbors and vice versa. Then one gets 
a subgraph of the tic-tac-toe graph such that each of the n 2 internal vertices lies 
on exactly 2 of the selected edges, and the An external vertices, taken in cyclical 
order, alternate between lying on a selected edge and not lying on a selected edge. 
Moreover, every such subgraph arises from an ice state in this way. 

Let us say that the leftmost vertex in the top row of external vertices is even. 
Figure 14 shows the seven subgraphs that result from applying the transformation 
to the seven ice states of order 3. 
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Figure 14: The seven FPL states of order 3. 

Leaving aside the behavior at the boundary, these are states of what physicists 
call the fully packed loop (FPL) model on the square grid (see e.g. [0]). I some- 
times prefer to call such states "near 2-factors" since nearly all of the vertices 
in these subgraphs have degree 2; only the external vertices have smaller degree 
(specifically, they alternate between having degree and having degree 1). 

If one starts from an external vertex, there is a unique path that one can follow 
using edges in the subgraph; this path must eventually lead to one of the other ex- 
ternal vertices. In addition to these paths ("open loops"), the edges of the subgraph 
can also form closed loops (see Figure 15, for example). 

Note that these loops (both open and closed) cannot cross one another. In 
particular, the open loops must join up the In external edges in some non-crossing 
fashion. If one numbers the vertices of degree 1 in cyclic order from 1 to In, 
the FPL state yields a pairing of odd-indexed external vertices with even-indexed 
external vertices. For instance, the FPL state shown in Figure 15 links 1 with 12, 
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2 with 11, 3 with 4, 5 with 6, 7 with 8, and 9 with 10. 




9 8 7 



Figure 15: A fully packed loop state of order 6. 

It has been conjectured, on the strength of numerical evidence, that the number 
of ASMs of order n in which the open paths link 1 with 2, 3 with 4, . . . , and 2n — 1 
with 2n is exactly the total number of ASMs of order n — 1. 

We do not have a proof of this, but curiously, we have a proof of something 
else: that the number of FPL states of order n in which the open paths link 1 with 
2, 3 with 4, . . . , and In — 1 with In is equal to the number of FPL states of order 
n in which the open paths link 1 with In, 2 with 3, 4 with 5, . . . , and In — 2 with 
2n — 1 . Note that when n is divisible by 4, the geometries of the two linking- 
patterns is different, with respect to the tic-tac-toe graph. Yet the number of FPL 
states is the same. 

This is a special case of a far more general fact proved by Wieland pB|]. For 
any two non-crossing pairings % and %' of the numbers 1 through 2n (viewed 
as equally spaced points on a circle), if it and it' are conjugate via a rotation or 
reflection, then (if we now treat the numbers 1 through 2n as the labels of vertices 
of degree 1 in the tic-tac-toe graph of order n) the number of FPL states with 
linking -pattern % equals the number of FPL states with linking-pattern %' . It is as 
if the tic-tac-toe graph, in some mystical sense, had an automorphism sending 1 
to 2, 2 to 3, etc. 

One might also ask for the number of FPL states of order n in which 1 is linked 
with 2 (ignoring all the other linking going on). Here too we have a conjectural 
answer, due to David Wilson: the number of such FPL states is just the total 
number of FPL states of order n multiplied by 

3n 2 + l 
2 4n 2 -l' 
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This would imply, in particular, that as n goes to infinity, the probability that a 
randomly chosen FPL state links 1 with 2 is asymptotically 3/8. 

For a very recent discussion of the FPL model, see pT|]. In this article, Razu- 
mov and Stroganov point out that the FPL model is closely related to a seemingly 
quite different lattice model. 

It is worth pointing out that these conjectures are truly native to the FPL in- 
carnation of ASMs; it is hard to see how they could have arisen from one of the 
other models. So, even though the transformation between ASMs and FPL states 
is fairly shallow mathematically, some interesting questions can arise from it that 
might not otherwise have been noticed. Likewise, the passage between square-ice 
states and ASMs is not conceptually deep, but it made possible the shortest known 
solution of the ASM problem by putting the problem into a form where known 
methods from the physics literature could be applied. Hence a good subtitle for 
this paper might have been "The non-trivial power of trivial transformations". 



8 Numerology 

As far as I know, the first manifestation of the sequence 1 ,2,7,42,429,7436,. . . occurred 
in connection with combinatorial objects called descending plane partitions or 
DPPs [|l]] . Another manifestation was totally symmetric self-complementary plane- 
partitions (TSSCPPs) Andrews discovered the proofs of both these formulas. 
Indeed, it was Andrews' proof of the 1 ,2,7,42,. . . formula for TSSCPPs that galva- 
nized Zeilberger into tackling the ASM conjecture. Zeilberger showed that ASMs 
are equinumerous with TSSCPPs. However, his proof was not bijective, and to 
this day nobody knows of a good bijection between ASMs of order n and TSS- 
CPPs of order n. 

Another context in which these numbers arise is the study of the XXZ model 
in statistical mechanics [fZDp [0- 

A different context in which numbers related to ASMs have occurred is certain 
"number walls" investigated by Somos p3p. (Number walls are arrays of Hankel 



determinants, arranged so as to faciliate calculations of successively larger ones; 
see []9|] for details.) In Somos' examples it is not the sequence 1,2,7,42,. . .that 
crops up, but sequences that enumerate various sorts of symmetric ASMs. Xin [J25p 
has found a Hankel determinant theorem that involves the sequence 1,2,7,42, .. . 
itself: he has shown that for all n, the determinant of the n-by-n matrix whose 
z,jth entry is equal to the coefficient of x' + i~ 2 in the Taylor expansion of the 
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generalized Catalan generating function — v 3x ' — [ |17| ] is equal to 

3© 

times the number of n-by-n alternating-sign matrices. 

It would be desirable to have some sort of understanding of why the num- 
ber of ASMs of order n (with or without symmetry-constraints) turns up in these 
seemingly disparate situations. 

9 Large random ASMs 

Another sort of phenomenon associated with ASMs of order n is their typical 
"shape" when n is large. I remarked above that the ASMs of order n form a 
distributive lattice; consequently, the method of "coupling from the past" can be 
applied [ |T9| ] for the purpose of generating a random ASM of fixed order governed 
by the uniform distribution. Figure 16 shows a random ASM of order 40, repre- 
sented as a gaskets-and-baskets tiling (where an attempt has been made to give 
each of the six tile-types its own distinctive shading). 

Note that the gaskets (which correspond to the non-zero entries of the ASM) 
stay away from the corners. Computer experiments strongly indicate that this is 
typical behavior: the probability of finding a non-zero entry close to one of the 
corners appears to be quite small. Another way of expressing this is in terms of the 
entries h i • of the height-function matrix. Say that a location (z, j) in a particular 
height-function matrix of order n is frozen if the height there is equal to either the 
maximum possible height that any height-function matrix of order n can exhibit 
at that location or the minimum possible height that any height-function matrix of 
order n can exhibit at that location. Then the claim is that a significant portion of 
the height-function matrix, concentrated in the four corners, tends to be frozen. 

This is analogous to phenomena that have been observed for other sorts of 
combinatorial models. Indeed, if one adopt a non-uniform distribution on the set 
of ASMs of order n, where the probability associated with an ASM containing 
exactly k entries equal to — 1 is proportional to 2 , then it is rigorously known that 
the frozen region tend in probability to a perfect circular disk [ [13| ] [|8|]. However, 
it is not known that this holds when the uniform distribution on ASMs is used. 
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Figure 16: A random gaskets-and-baskets tiling of order 40. 
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10 Back to Dodgson 



I conclude this article by coming full circle and returning the context in which 
ASMs first came to light: the study of Dodgson's condensation algorithm and 
its variants. I will not discuss Dodgson's algorithm per se, but rather a variation 
of it invented by Robbins and Rumsey []24|]. This modified form of Dodgson 
condensation is an algebraic recurrence relation 

fi,j,k+l = (fi+lj,kfi-l,j,k + fi,j+l,kfij-l,k)/ fi,j,k-l 

satisfied by certain functions / : Z 3 — > R. (Note that this equation can be written 
slightly more symmetrically as 

fi,j,k+lfi,j,k-l ~ fi+lj,kfi-lj,k~ Aj+ l,kAj-l,k = ®- 

Physicists and researchers in the field of integrable systems call this a discrete 
Hirota equation, and have developed a great deal of theory associated with it; 
however, the observations I make here seem to be currently unknown outside of a 
small circle of algebraic combinatorialists.) If we let /. . _j = x i ■ and f i j Q = y t ■ 

for formal indeterminates x i ■,y i ■ (with (z, j) ranging over Z 2 ), then the recurrence 
relation lets us express all the /• . ^'s in terms of the x t and y i y's, at least for- 
mally. A priori, one expects each /. . k to be a rational function of the x and y 
variables; the surprise (the first of several surprises, in fact) is that these rational 
functions are actually Laurent polynomials (that is, they are polynomials functions 
of the x and y variables along with their reciprocals). 

This observation seems to have first been made by Mills, Robbins and Rumsey, 
who actually considered a more general recurrence 

fij,k+l = (fi+l,j,kfi-l,j,k + ^fi,j+l,kfi,j-l,k)/fi,jJc-V 

The case A = — 1 corresponds to the original Dodgson algorithm, but Mills et al. 
noticed that the same surprising cancellation occurs for more general values of 
A, including the especially nice case A = 1, and that one always obtains Laurent 
polynomials. 

It was from studying these Laurent polynomials that Mills, Robbins and Rum- 
sey were led to discover alternating- sign matrices. Every term in one of these 
Laurent polynomials has a coefficient equal to 1 (that is the second surprise), and 
is a product of powers of a finite number of x and y variables. The third surprise 
is that all the exponents of the variables are +1, —1, and 0. The fourth and final 



17 



surprise is that these patterns of exponents encode ASMs. More specifically, the 
exponents of the x variables (after a global sign flip) encode one ASM, and the 
exponents of the y variables encode another. These two ASMs satisfy a combina- 
torial relationship that the researchers dubbed "compatibility". They showed that 
the number of compatible pairs of ASMs is exactly 2"(" +1 )/ 2 . 

As it happens, this formula is not hard to verify as a consequence of the other 
claims I have made. If indeed all the coefficients in the Laurent polynomial equal 
1, then one can count the terms (and thereby count the compatible pairs of ASMs) 
just by setting all the x and y variables equal to 1 . But in this case, /• j k depends 
only on k (call it F,), and the three-dimensional recurrence boils down to the one- 
dimensional recurrence 



F k + l = ( F k F k + F k F k)/ F k-l 

with initial conditions F Q = F x = 1, which is readily solved. 

The terms of these Laurent polynomials were originally understood in terms of 
compatible pairs of ASMs. A few years after this work was done, it turned out that 
compatible pairs of ASMs admit a much more geometrical representation, namely, 
as tilings of regions called Aztec diamonds by means of tiles called dominos (1- 



by-2 and 2-by-l rectangles). See Ql 1Q for more details. 

I will close by pointing out that a kindred recurrence relation cries out to be 
studied, namely 

fi,j,k = (fi-lj,kfij-l,k-l +fi,j-Lkfi-l,j,k-l +fij,k-lfi-l,j-l,k)/ fi-Lj-Lk-V 
with initial conditions 

[ x ij, k if = -i, 

[ 'i.j.k ifi + j + k = +l. 

Here (just as in the Mills, Robbins, and Rumsey recurrence) one finds empirically 
that each value of /. . k is expressible as a Laurent polynomial in the x, y, and z 
variables (in fact, shortly before this article went to press, this Laurent property 
was proved by Fomin and Zelevinsky[JT2p); here too one finds empirically that 
each coefficient in these Laurent polynomials is equal to 1 ; and here too one finds 
that the exponents of the x, y and z variables that occur in the Laurent monomials 
are universally bounded (in this case between — 1 and +4 rather than between — 1 
and +1). If all this is true, then the exponent-patterns that arise are some sort 
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of analogue of compatible pairs of ASMs, and moreover, we know exactly how 
many there are: 3 L n / 4 J . (This comes from reducing the original three-dimensional 
recurrence to a one-dimensional recurrence, as we did before.) So, assuming that 
our empirical observations are not leading us astray, there is some new kind of 
combinatorial gadget that governs these Laurent polynomials (or vice versa!), and 
we know exactly how many gadgets of order n there are. And it is easy to generate 
these Laurent polynomials (and with them the gadgets) using Maple, e.g. with 
the following short program: 

f : = proc (i, j, k) 

if (i+j+k < 3) then x(i,j,k) else 
simplify ( 

( f (i-1, j,k) *f (i, j-l,k-l) + 

f (i, j-l,k) *f (i-1, j,k-l) + 

f (i, j,k-l) *f (i-1, j-l,k) ) 
/f (i-1, j-l,k-l) ) ; 
fi; end; 

Nonetheless, we do not know what these gadgets are, combinatorially! They are 
analogous to pairs of compatible ASMs, which in turn are equivalent to domino 
tilings of Aztec diamonds, so one hopes that the gadgets have some intuitive geo- 
metric meaning. For more information about the properties of these gadgets, see 
http : / / www . math . wise . edu/ ^propp/ cube-recur . 

More late-breaking news: Harvard undergraduates Gabriel Carroll and David 
Speyer solved this problem in 2002, as participants in the Research Experiences 
in Algebraic Combinatorics at Harvard program. More details will at some point 
become available through links at the REACH web-site: 

http : / / www . math . harvard . edu/ ^propp/ reach/. 
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